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Abstract We calculated the glassy response of solid 4 He to torsional oscillations assuming 
a phenomenological glass model. Making only a few assumptions about the distribution of 
glassy relaxation times in a small subsystem of otherwise rigid solid 4 He, we can account 
for the magnitude of the observed period shift and concomitant dissipation peak in several 
torsion oscillator experiments. The implications of the glass model for solid 4 He are three- 
fold: (1) The dynamics of solid 4 He is governed by glassy relaxation processes. (2) The 
distribution of relaxation times varies significantly between different torsion oscillator ex- 
periments. (3) The mechanical response of a torsion oscillator does not require a supersolid 
component to account for the observed anomaly at low temperatures, though we cannot rule 
out its existence. 
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1 Introduction 

Torsion oscillators have been used successfully to measure an anomalous change in resonant 
period and accompanying dissipation in solid 4 He. 'i 2 ! 3 ! 4 ! 5 ! 6 ? 7 It has been speculated that the 
observed signature is due to Bose-Einstein condensation of vacancies or interstitials form- 
ing a novel supersolid state in otherwise crystalline 4 He. 8 i 9 i 10 i 11 i 12 i 13 Early on, the change 
in resonant period has been attributed to a (nonclassical) decoupling of a supersolid compo- 
nent from the (classical) normal moment of inertia. This is not surprising, since the observed 
change in period is in agreement with similar observations of onset of superfluidity in liquid 
4 He, measured a long time ago by torsion oscillators. 14 i 15 i 16 i 17 However, it is important to 
remember that for liquid 4 He it was already well established, long before the torsion oscilla- 
tor (TO) experiments were performed, that it undergoes a transition from liquid to superfluid 
with no viscosity. It was natural to use the connection between superfiow and period drop. 
Status of a search of supersolidity in solid 4 He is different: a change in period has been re- 
ported, but no evidence of mass superfio w 18 - 1 9-20,21 1 22 1 23 1 24 Qr condensation— has been 
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Fig. 1 (Color online) Hysteresis of rim velocity of nonclassical rotational inertia fraction (NCRIf) by Aoki et 
al.— Measurements were taken at 19 mK as the oscillator drive was increased (circle) and decreased (square). 



seen below the expected transition temperature. It is therefore necessary to ask the ques- 
tion what is the relationship between the change in period and superflow, and what alternate 
physical mechanisms can explain the change in period and concomitant peak in dissipation. 

A potential contender for an alternate explanation is the glass scenario, which automati- 
cally accounts for hysteresis effects, annealing dependence, linear term in specific heat, and 
long relaxation times in many observables. For example, the strong hysteresis effect reported 
by Aoki et al.— for the rim velocity dependence of the torsion oscillator, shown in Fig. [T] 
is consistent with a glassy response. However, it is difficult to explain within a supersolid 
scenario how increasing the rim velocity of the cell does not change the reported nonclas- 
sical rotational inertia: between slow and fast rim velocity one increases kinetic energy of 
solid, °c v 2 , and the ratio of kinetic energies is on the order of ~ (600/10) 2 = 3600. Fast 
rim velocity exceeds the supposed critical velocity of the condensate, ~ lOjUm/s, by several 
orders of magnitude. 

In this article, we explore if a phenomenological glass model can account for the change 
in period without postulating a nonclassical moment of inertia. Similar to the supersolid pic- 
ture, we assume that only a small subsystem of solid 4 He exhibits glassy properties that dom- 
inate the response at low temperatures. This is an important point, since it has been argued 
before that the observed large change in dissipation cannot be described by uniform Bose- 
Einstein condensation. 28 i 29 i 30 It remains to be seen if nonuniform Bose-Einstein condensa- 
tion alone either along grain boundaries 31 or along the axis of screw dislocation s 32 i 33 i 34 can 
explain the dynamic response of TOs. A discussion on the role of a glassy component does 
not rule in or out the presence of a supersolid component. We emphasize that our analysis 
addresses a glassy contribution regardless of the magnitude of a supersolid component. 



2 Glass model for torsion oscillator 



Glass model: In previous work, we argued for the possibility of a glass phase at low tem- 
peratures, roughly below ~ 150 mK, to explain the observed anomalous linear temperature 
dependence in the specific heat of the otherwise perfect Debye solid 4 He. 35 ' 29 Below we will 
give an extended oscillator analysis for which the exact nature of the glass is not crucial. For 
example, it may be caused by two-level systems (TLS) of pinned dislocation lines, vortex ex- 
citations, etc. However, it is important to point out that the amplitude of the period shift and 
dissipation peak can be changed dramatically by the growth history or annealing process of 
the crystal. In order to explain these puzzling features of solid 4 He, we conjectured earlier— 
that structural defects, e.g., localized dislocation segments, form a set of TLS observable 
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at low temperatures. These immobile crystal defects affect the thermodynamics 35 ' 9 of bulk 
4 He and the mechanics 2 * of the TO loaded with 4 He. For the analysis of the specific heat, 
we used independent TLS to obtain the universal glass signature of a linear-in-temperature 
specific heat term at low temperatures. In parallel, we used a phenomenological glass model, 
that may originate from an ensemble of TLS, to describe the mechanical response of the TO. 

Torsion oscillator and rotational susceptibility: To set the discussion, we note that TO 
experiments measure the period and dissipation. One applies a force and generates a dis- 
placement of the oscillator. The relationship between the force and displacement (angle) is 
controlled by the TO susceptibility. In order to extract information from such an experiment, 
one needs a model to determine the relation between observable period and dissipation of 
the TO and the corresponding moment of inertia, damping and effective stiffness of the me- 
dia. We start with the general equation of motion for a harmonic TO defined by an angular 
coordinate Q in the presence of an external and internal torque,— 

Io S cO(t) + Y osc O(t) + a osc e(t) =M ext (t)+M int (t). (1) 

Here, I osc is the moment of inertia of the (empty) TO chassis, a osc is the restoring (stiffness) 
coefficient of the torsion rod, and y osc is its dissipative coefficient. M ext (t) is the externally 
imposed torque by the drive. M- mt (t) = J g(t — t')d(t')dt' is the internal torque exerted by 
solid 4 He on the oscillator for a system with time translation invariance. In general, the 
backaction g(t — t') is temperature, T, dependent. The experimentally measured quantity is 
the angular motion of the TO - not that of bulk helium, which is enclosed in it. Ab initio, 
we cannot assume that the medium moves as one rigid body. If the non-solid subsystem 
"freezes" into a glass, the medium will move with greater uniformity and speed. This leads 
to an effect similar to that of the nonclassical rotational moment of inertia, although its 
physical origin is completely different. Therefore, we argue for an alternate physical picture, 
namely that of softening of the oscillator's stiffness. The angular coordinate d(t) of the TO 
is a convolution of the applied external torque M ext (t) with the TO susceptibility %(t,t'). 
Under Fourier transformation the angular response of the TO is 

X Q l {co)e{co)=M exl {co) +M ini (w). (2) 

Defining the total angular susceptibility as %~ l = Xq 1 ~ Mj„ t , we write 

X~\g>) = a osc -iy osc co-I osc co 2 -g{co), (3) 

where g(co) is the Fourier transform of the backaction due to the added solid 4 He. In what 
follows, we will treat the backaction as a small perturbation to the TO chassis. 



3 Period and dissipation of torsion oscillator 

We now determine the experimental consequences of the phenomenological glass model, 
where a small glassy subsystem of solid 4 He gives rise to the observed dynamic behavior. 
Glass is generally defined as a supercooled liquid out of equilibrium on measurable time 
scales: its equilibration time becomes extremely large (and unmeasurable) at low temper- 
atures. Any glass former is a liquid at high temperature and becomes an amorphous solid 
(the glass) at sufficiently low temperatures. In this context bulk solid 4 He is not a glass for- 
mer; we are talking about glass forming within a small fraction of 4 He sample. We note that 
our analysis will remain qualitatively unchanged for a general description of the system by a 
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"freezing" at low temperature of the appropriate component (defect or other) that is dynamic 
at high temperatures (see the appendix of Nussinov et al.— ). 

We start by reviewing results for an underdamped harmonic torsion oscillator. The reso- 
nant period is obtained from the angular coordinate 9(t) = Re{0oexp[— icot]}, with a com- 
plex amplitude 9q and complex frequency CO = COq — iK. In the case of an underdamped TO 
with K <C COq, the resonant period is P = 2k/<j}q, and the quality factor Q or dissipation is 
Q _i = 2k/ COq, with resonant frequency COq = \J a osc /I osc . 

In the remainder, we use effective oscillator parameters, which are defined as the sum of 
parameters describing the chassis, x^ 1 , and the added solid 4 He given by 

g(ffl) =iy H eCO+lHeC0 2 +&((»). (4) 

Thus, we write the net moment of inertia I = I osc + In e and dissipation y = y osc + yu e . The 
total response function of the TO is given by Eqns. (13} and @J, The term &(co) captures 
the dynamics of a glass component and is a function of temperature and frequency. In the 
limit co — > the term &(co) — > as the mechanical motion of any glass component will 
be the same as that of the surrounding solid. Hence there will be no relative motion and 
no transient overdamped modes for co = 0. However, at any finite frequency co, we can ap- 
proximate the glass response by ( S{co) « goG(co), where the coefficient go measures the 
glassy contribution of the solid and is evaluated at the resonant frequency coo of the TO. 
The dynamic response function G(co) of a glass can be approximated by a distribution of 
overdamped oscillators with different relaxation times T. Two popular relaxation time dis- 
tributions used in the literature are the Cole-Cole (CC) and Davidson-Cole (DC) functions. 
Both describe a superposition of overdamped oscillators. 36 ' 37 The CC distribution gives 
G(co) = 1/[1 - (icox)"], while the DC distribution results in G(co) = 1/[1 - icox]P . 

By comparison to Eq. ® for the TO chassis system with no helium, the glassy part of 
the backaction of 4 He, ^(co), renormalizes the effective spring stiffnes s 28,38 

a eff ~ (a osc - go) , for cot < 1 , (5) 
a eff ~ a osc , for cot > 1. (6) 

These expressions flesh out the dependence of the medium response on the applied driving 
frequency. When the driving frequency is far more rapid, co ^> T _1 , then the transient re- 
sponse of the medium is that of a liquid. In that limit, the transient modes within the medium 
cannot "keep up" with the driving torque and only the bare stiffness of the TO remains aug- 
mented by the solid helium contribution. The effective spring stiffness is that of the driving 
oscillator, a e ^ = CC osc , see Eq. (5}. The limit T~' — > corresponds to that of an ideal rigid 
low-temperature glass in which no transient liquid-like response of the system is present. 
By contrast, for slow oscillations co <?C T~', the excited modes within 4 He are of charac- 
teristic transient time T that is long enough to respond to the driving torque and lead to an 
additional backaction and effective reduction of the spring stiffness, see Eq. ©. From this 
discussion it is clear that the maximum relative shift in period or frequency will depend on 
the glassy fraction go given by Aco max /cOo ~ go/(2a OJC ), which can vary widely between 
different torsion oscillators, growth and annealing procedures. 

The resonant frequency of the TO with backaction is given by the root of 

X~\co) = a-iyco-Ico 2 -g G(co) = 0. (7) 

We anticipate that when the relaxation time is similar to the period of the underdamped TO, 
the dissipation will be maximal. Here, the glassy component responds with the same fre- 
quency as the "normal" solid component. The glassy part merely renormalizes the effective 
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spring constant a, but does not lead to additional transient modes, which closely interfere 
with the oscillations of the "normal" part of the TO. We look for the largest magnitude of the 
imaginary part of the root and see when it is maximal as a function of T. A larger imaginary 
part implies a shorter decay time and a smaller value of Q . Since the homogeneous Eq. I0 
is scale invariant, we normalize all oscillator quantities by the effective moment of inertia /, 
i.e., a = a //, y = 7//, and | = go/I- 

As can be seen from Eq. Q, for an ideal dissipationless oscillator, 7 = 0, the resonant 
frequency = \/~G. is the pole of %{(0) in the limit T — > 0. If we expand X~ l about this 
root, co = a>o + 8co, with 8 CO = C0\ — iK, then we find to leading order in 8 CO 

8co^- i?COo + ^ G{COo) . (8) 
iy+2coQ 

Therefore, the root attains an imaginary component K and the dissipation becomes 

e _1 = — « — ImG^ + g- 1 , (9) 

COq COo 

with A = go/coo and Q~ l = y/coq. As COq increases for fixed a osc , Q^} increases. For a ~ 
B ~ 1, the dissipation peaks near coqT ~ 1. Similarly the resonant frequency becomes 

2k A / \ 

2k f= — ^co -—{2cooReG{coo)+flmG{coo)j, (10) 

which increases monotonically when T is lowered. For the special case of Debye relax- 
ation processes ( a = /3 = 1), we find Re G(coq) = [1 + (<Bot) 2 ] -1 and Im G(coq) = C0qt[1 + 
(cOqT) 2 ] -1 and recover results reported earlier,— except for the additional contribution pro- 
portional to 7 in Eq. d 1 Ob - It follows that the changes in dissipation, AQ~ l = Q~ l — Q~\ and 
frequency, Aco = COq — 2n/P, determine the glass relaxation time T. Combining Eqns. ^ 
and < l 1 Ob we arrive at a general relation between shift in dissipation and frequency for fx <C 1 : 



AQ- 1 4ImG 2 Im G 

Aco ~ 2ct)oRe G + flm G ~ coo Re G' 

For example, for a DC glass distribution this becomes 
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AO' 1 2 f 2/3T ' G * )T " ' 

~ _ tan (B arctan(ffloT)) ~ <^ tan(jS^/4), cOqT = 1, (12) 
Aco coo ™, ,n ,„s 

t ^tan(jS^/2), coqT^oo. 

In the past, there have been several reports of large experimental ratios A Q 1 ^ ~ 3 — 
j 2 2^4,5 B ecause of jS < 1 and for COqT ~ 1 the ratio is limited to AQT X < 2, this requires 
diverging relaxation times close to the temperature where the dissipation peaks. For such 
cases, COqT — > °° in Eq. d 1 2b and ratios of order 10 can be obtained for values of B ~ 0.6. 
On the other hand, for B = 1, Eq. i fTTt simplifies even further with AQT 1 x2cooi:(Aco/coo). 
Similar results for the ratio were obtained for other phenomenological models with dissipa- 
tive channels. 30 ' 38 For example, Huse and Khandker 30 assumed a simple phenomenological 
two-fluid model, where the supersolid is dissipatively coupled to a normal solid resulting in 
a ratio of AQ^ 1 ^ « 1. Yoo and Dorsey 38 developed a viscoelastic model and Korshunov— 
derived a TLS glass model for solid 4 He that captures the results of the general phenomeno- 
logical glass model originally proposed by Nussinov et al.— Here, we like to emphasize 
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Fig. 2 (Color online) The period shift AP = P — Pq (black, left axis) and dissipation (red, right axis) vs. 
temperature for solid 4 He. The experimental data are from Rittner and Reppy, Fig. 2 of Ref A A Davidson- 
Cole (DC) and Cole-Cole (CC) fit are shown. The DC fit was performed with parameters j3 = 0.60, = 
1 1.4- 10~ 6 , f = 184.2305 Hz, 5/ = 69 /iHz, A = 33.8 mHz, T = 0.439 /is, DT = 1.173 K, T = K, a T = 
2.0- 10~ 5 KT 1 . The CC fit used parameters a = 1.15, QZ 1 = 11.1 ■ 10~ 6 , fo = 184.2305 Hz, Sf = 20 /iHz, 
A = 34.7 mHz, T = 1.95 tis, DT = 0.868 K, T = K, a T = 1.6- HT 5 K _1 . 

that it is challenging to reconcile a large dissipative AQ/Q ratio with uniform Bose-Einstein 
condensation.— ' 29 ' 30 

We now make further assumptions about the glassy relaxation time T. In many glass 
formers T follows the phenomenological Vogel-Fulcher-Tamman (VFT) expression t{T) = 
Toexp[D^o/(7 , — To)] for T > Tq. Here, To is the temperature at which an ideal glass tran- 
sition occurs, which is below the temperature where the peak in dissipation occurs. The 
parameter D is a measure of the fragility of the glass (D < 10 for fragile glasse s 40 ' 41 ). Fi- 
nally, at temperatures T < To the glassy subsystem freezes out and T becomes infinite. 



4 Results and discussion 



All samples of solid He used in this study had in common that they were grown with the 
blocked capillary method using commercial grade helium ( 3 He impurity level ~ 0.3 ppm). 
Also, it is important to remember that both glass models (CC and DC) use only five fit pa- 
rameters: go, To, DTq, To, and either an exponent a or j8. All other oscillator parameters 
are determined by normal state properties of the TO loaded with solid 4 He. In addition, 
we noticed during our analysis of the TO experiments that in order to fit the glass mod- 
els to the experimental data, we had to correct the resonant frequency by a small amount, 
f = fo + Sf, because in many reports fo is not available to desired absolute accuracy or 
data are only reported relative to a high-temperature resonant frequency. Furthermore, the 
fits were complicated by the experimental observation of a slight temperature dependence 
of the resonant frequency at higher temperatures. To account for this drift in frequency of 
yet unknown origin, we approximated fo(T) « /q(0)[1 — UjT] by a small linear-in-T cor- 
rection. 
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Fig. 3 (Color online) The resonant frequency (black, left axis) and dissipation (red, right axis) vs. tem- 
perature. The experimental data are from Hunt et al. 7 with Cole-Cole (CC) parameters a = 1.85, Q~ l = 
1.23 ■ 1(T 6 , f = 574.4768 Hz, Sf = 1.489 mHz, A = 347 mHz, T = 2.52 /is, DT = 0.408 K, 7b = -44 
mK, a T =2.43- 10~ 5 K _1 . 



The TO experiment reported by Rittner and Reppy, 2 see Fig. [2] is in excellent agreement 
with the proposed glass models. Both CC and DC glass distributions require exponents 
different from unity, which means that there is a spread of relaxation times T. 

In Fig. [3j we report an analysis of the measured data by Hunt et al. 7 assuming a CC 
distribution of relaxation times. As can be seen, an excellent fit is obtained. For comparison, 
we also tried a DC distribution for relaxation times, but found only fair agreement. It is worth 
pointing out that unlike in the Debye relaxation analysis by Hunt and coworkers (a single 
overdamped mode), we do not require a supersolid component to simultaneously account 
for frequency shift and concomitant dissipation peak. 

Finally, in Fig. |4] we report a DC and CC analysis of the measured data by Aoki et 
al.— for the in-phase mode of their double resonance compound TO. We obtain excellent 
agreement between experiment and glass model assuming a CC distribution, while a DC 
distribution for glassy relaxation times results only in fair agreement. 



5 Conclusions 

To summarize, we have shown that a phenomenological glass model describing a small 
subsystem of solid 4 He can simultaneously account for the experimentally observed change 
in resonant period (frequency) and the concomitant peak in dissipation. 

Our analysis of TO experiments reveals that most are better described by a Cole-Cole 
distribution for glassy relaxation times. Unlike for conventional structural or dielectric glasses, 
where the CC exponent a is usually less than unity, we find consistently a > 1. This may 
reflect on the possible nature of a quantum or superglass in solid helium. Further, we derived 
a simple relation for the ratio of change in dissipation and change in resonant frequency (pe- 
riod) that can explain the large ratios of order ~ 10 observed in experiments. The values for 
the glass exponents a or jS required to fit the experiments by the Rutgers and Cornell groups 
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Fig. 4 (Color online) The period shift (black, left axis) and dissipation (red, right axis) vs. temperature. The 
experimental data are from the in-phase mode (495.8 Hz) of the coupled double oscillator by Aoki et al^ 
The experimental data are already corrected for temperature dependence. The DC parameters are j8 =0.12, 
Q- 1 = 0.824 ■ 10- 6 , f = 495.829 Hz, Sf = 4.4 /iHz, A = 8.19 mHz, T = 2. 15 /is, DT = 0.306 K, 7b = 23.6 
mK, ot r =0K-'.TheCC parameters are a = 1.70, Q' 1 = 0.793 ■ 10~ 6 , f = 495.829 Hz, 5/ = -19 /iHz, 
A = 8.97 mHz, T = 13.2 /Is, £>7/o = 0.248 K, T = -17 mK, a r = KT 1 . 

point toward broad distributions of glassy relaxation times. This invalidates any attempt to 
describe these experiments by a single overdamped mode (Debye relaxation). These glassy 
relaxation processes should also have significant effects on thermodynamics and dynamics 
of solid 4 He. The key result of this work is that many TO experiments can be described 
assuming that a small fraction of solid 4 He undergoes a glass transition at low temperatures. 
Whether or not there is a supersolid fraction present in solid 4 He is beyond this analy- 
sis. A frequency-tunable TO may differentiate between a glassy contribution leading to an 
increase in the maximum frequency shift, Aco max ~ (OQgo/oc osc , and no change in the dissi- 
pation shift, zlg 1 ~ go/oiosc, with increasing Oo, while the frequency shift for a supersolid 
should decrease with increasing <x>q. Our study shows that the unequivocal identification of 
supersolidity in solid 4 He is challenging and does require clear understanding of normal state 
dynamics. Clearly, more dynamic studies probing the frequency or time response to a stimu- 
lus and detailed bulk characterization of samples are necessary to investigate the differences 
between small subsystems of glassy, supersolid or superglassy origin. 
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